Excitation Spectrum and Superexchange Pathways in the Spin Dimer 

VODPO4 I/2D2O 
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Magnetic excitations have been investigated in the spin dimer material VODPO4 ■ I/2D2O us- 
ing inelastic neutron scattering. A dispersionless magnetic mode was observed at an energy of 
7.81(4) meV. The wavevector dependence of the scattering intensity from this mode is consistent 
with the excitation of isolated V*"^ spin dimers with a V-V separation of 4.43(7) A. This result is 
unexpected since the V-V pair previously thought to constitute the magnetic dimer has a separation 
of 3.09 A. We identify an alternative V-V pair as the likely magnetic dimer, which involves superex- 
change pathways through a covalently bonded PO4 group. This surprising result casts doubt on the 
interpretation of (VO)2P207 as a spin ladder. 



PACS numbers; 75.10.Jm, 75.30.Ds, 75.30.Et 

One of the goals of modern condensed matter physics 
is the comprehensive understanding of the ground state 
structure and excitations in quantum magnetic systems. 
Low dimensional systems, including spin chains, planes, 
and intermediate dimensional spin ladders are of partic- 
ular importance because of the role they may play in the 
physics of high-Tc superconductors The exchange 
coupled spin dimer is the simplest interacting quantum 
magnetic system. Spin dimers can be thought of as 
the basic building blocks for more complicated spin sys- 
tems, including RVB spin liquids and spin ladders. In 
this Letter we present the results of a neutron scatter- 
ing study of excitations in a dimer system of particu- 
lar interest: deuterated vanadyl phosphate demihydrate, 
VODPO4 • I/2D2O , (abbreviated VODPO). VODPO is 
closely related to the compound (VO)2P207 (VOPO) 
which is widely thought to be an excellent reali- 
sation of the two-leg antiferromagnetic Heisenberg spin 
ladder. The present results show that there has been a 
fundamental misunderstanding of the magnetism in these 
materials, and that a careful re-examination of conclu- 
sions drawn about VOPO is called for. In particular the 
dominant exchange couplings have been misidentified. 

The spin ladder nature of VOPO was proposed be- 
cause the crystal structure contains S — 1/2 V^"*" ion 
"rung" pairs which are stacked into ladders Mag- 
netic interactions were thought to occur through V-O-V 
superexchange paths along both rung and chain direc- 
tions. Neighbouring ladders are separated by PO4 com- 
plexes which have been assumed to be magnetically in- 
ert. (VO)2P207 can be synthesized by dehydration of the 
precursor material VOHPO4 • I/2H2O (VOHPO) |,|. 
VOHPO contains V-V pairs identical to the proposed 



rungs of the VOPO ladder, but magnetically decoupled 
along the chain direction (see Fig. 1). It is therefore pos- 
sible to study part of the ladder in isolation in VOHPO 
or the deuterated counterpart VODPO. For neutron scat- 
tering experiments the deuterated material has the ad- 
vantage of lowering the incoherent background without 
significantly altering the magnetic interactions. 

An idealized isotropic spin-dimer material consists of 
isolated pairs of magnetic ions that interact through a 
Heisenberg Hamiltonian, 



H — J Si ■ S2 
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Johnson et al. measured the magnetic susceptibility 
of VOHPO and found that it is accurately described by 
the Bleaney-Bowers formula for the susceptibility of iso- 
lated dimers of = 1/2 spins, with J = 7.6 meV and 
g = 1.99 (close to the free electron value, as expected 
for an orbitally quenched V^+ ion). 

By measuring the E and Q dependence of magnetic 
excitations using inelastic neutron scattering we can test 
this simple Hamiltonian and deduce the distance between 
V-V ion pairs in a dimer. Since a dimer has an S total = 
ground state and only one excited level (an Stotai = 1 
triplet at Ei ^ Eq — J), inelastic neutron scattering 
should reveal a single, dispersionless excitation. The sus- 
ceptibility fit predicts that this excitation should appear 
at Egap = 7.6 meV. 

On carrying out a powder average over the neutron 
scattering structure factor S{Q, E) of a dimer, the pow- 
der scattering intensity I{Q, E) is found to be 



I{Q,E)^\F{Q)\H^- 



sin QR 
QR 



5{E - Egap) . (2) 
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Although we expect only a single excited level, the Q de- 
pendence of the magnetic excitation is especially inter- 
esting in that it allows a determination of the separation 
R of the magnetic ions forming the dimer pair. 

The VODPO used in this study was synthesised from 
V2O5, D3PO4, and D2O, which were refluxed under ni- 
trogen for 2 days to obtain VOP04-2D20. This prod- 
uct was heated in a Teflon-lined reaction vessel, with 
isopropanol-d8 under nitrogen and autogenous pressure 
at 125° for 7 days, and then was filtered and dried in 
a nitrogen glove bag to obtain VODPO. Infrared spec- 
troscopy indicated an H/D ratio of 0.04. 

Neutron powder diffraction was used to refine the crys- 
tal structure, which gave parameters similar to previous 
results for VOHPO Rietveld refinement yielded an 
orthorhombic structure at T = 10 K, with lattice pa- 
rameters of a = 7.4102(6) A, = 9.5861(8) A, and 
c — 5.6873(7) A. It also confirmed the atomic fraction of 
deuterium as being at least 96%. The crystal structure of 
VODPO at 10 K is illustrated in Figure la. The structure 
contains V*+ ions, each of which has a single 3d electron 
and effective spin S — 1/2. Neighboring V^"*" ions belong 
to edge-sharing VOe octahedra. Next-nearest-neighbor 
V^^ ions are separated by PO4 tetrahedra. 

Possible exchange paths between neighboring V*"*" ions 
are shown in Figure lb. The nearest-neighbor pathway 
(PI) has a V-V separation of dyv — 3.09 A and the 
next-nearest-neighbor pathway (P2) has a V-V separa- 
tion of dyv — 4.32 A. Consideration of orbital overlaps 
led Beltran-Porter et al. to conclude that superex- 
change in VODPO is probably dominated by the PI 
pathway. 

The magnetic excitation spectrum of VODPO was 
measured using the HBl triple- axis spectrometer at the 
High Flux Isotope Reactor at Oak Ridge National Labo- 
ratory. Approximately 10 grams of VODPO powder was 
sealed in an aluminum cylinder and mounted in a dis- 
plex cryostat, allowing the temperature to be varied over 
the range 10 K — 300 K, with an accuracy of ±0.1 K. 
The Q resolution of the spectrometer was relaxed be- 
cause the excitations were expected to be dispersionless. 
A vertically focussed PG(002) monochromator and hor- 
izontally focussed PG(002) analyser were used, and an 
instrumental coUimation of 48'-40'-open-open was cho- 
sen. Measurements were made using a fixed final energy 
Ep = 13.5 meV, giving an energy resolution of 1.16 meV 
full- width half maximum (FWHM). 

Scans in E at fixed Q were carried out in order to 
measure the magnetic excitations. Figure 2 shows the 
scattering at temperatures of 10 K, 50 K and 200 K, at 
a constant wavevector transfer oi Q = 1.0 A^^. In each 
scan there is a well-defined peak near 7.8 meV, close to 
the value of 7.6 meV anticipated from the susceptibil- 
ity data To extract accurate positions and widths 
for the peaks, a least-squares routine was used to fit a 
Gaussian profile with a linearly-sloping background. This 



gave a peak center of J = 7.86(6) meV and FWHM of 
2.05(10) meV for the scan at 10 K. Peak centers, intensi- 
ties and widths were extracted at the other temperatures 
using the same procedure. The integrated intensity falls 
with increasing temperature roughly as expected for a 
magnetic dimer mode. 

To investigate the dispersion of this mode we carried 
out a series of scans at 10 K with increasing wavevector, 
covering the range Q — 1.0-4.5 A~^, at intervals of AQ = 
0.1 A-^ Each scan was fitted to a Gaussian with linear 
background. The mode was found to be essentially non- 
dispersive, although the fitted width increases slightly 
with Q. 

To find the dimer spacing dyvj a constant energy scan 
was made at 7.75 meV. Figure 3 shows this scan at a 
temperature of 10 K. There are two clear maxima in the 
data, near 1.0 A~^ and 2.3 A~^; this oscillatory behavior 
is expected for scattering from a dimer, c/. Eq. (2). The 
solid line in the figure is a least-squares fit to Eq. (2), in 
which the (Jo) approximation for the form factor j|] 
has been used. A sloping background has been added to 
take account of the phonon contribution, which increases 
with Q. The fit yields a V-V separation in the dimer of 
dyv = 4.43(7) A, which we consider our definitive value. 

Figure 4 shows a two-dimensional color plot of the com- 
bined constant-Q scans at 10 K. The data have been 
interpolated to smooth the image, and the intensity is 
indicated by a colorbar. The spin excitation is observed 
as a band at 7.8 meV, with maximum intensity at the 
lower wavevectors. Two maxima of the structure factor 
as a function of Q appear as bright regions near 1.0 A~^ 
and 2.3 A^^. The non-magnetic phonon background is 
visible at large Q, especially at energies above 8 meV. A 
simultaneous least-squares fit was made to Eq. (2) plus a 
background term linear in E and Q. The delta function 
in energy was replaced by a Gaussian lineshape. The 
global fit gave an excellent description of the data with 
J = 7.81(4) meV, dyy = 4.49(8) A, consistent with the 
4.43(7) Afound in the constant-i? scan in Figure 3. A 
goodness-of-fit reduced = 1.12. The fitted model is 
shown as a contour plot (solid lines) overlying the data 
in Figure 4. 

Our constant- measurement gave a dimer spacing of 
dyy = 4.43(7) A, consistent with the exchange path P2 
shown in Figure lb (dyy — 4.32 A), but not with the 
shorter PI exchange path (dyy = 3.09 A) assumed to be 
dominant in earlier work . This is surprising because 
it implies that the exchange is stronger through a V-O-P- 
0-V pathway than through a much shorter V-O-V path- 
way. The conventional understanding of superexchange 
is that the interaction decreases rapidly with the number 
of ions in the exchange path |^ . We argue however that 
when there is a covalent complex present, such as the 
PO4 in VODPO, coherent orbitals can provide a strong 
exchange route. Understanding the electron wavefunc- 
tions within such covalent groups is therefore crucial in 
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predicting the important superexchange pathways. Pre- 
sumably the V-O-V superexchange is also weakened by 
the large departure of the bond angle (experimentally 
97.5°) from the optimum 180° for -0-. 

Beltran-Porter et al. have considered exchange 
pathways in several vanadyl phosphates. Although they 
did not anticipate a dominant PO4 pathway in VODPO, 
they do argue that PO4 pathways may be significant 
in (VO)2P207, as well as V-O-V. If V-O-V (rung) and 
V-O-P-O-V are the important exchange pathways in 
(VO)2P207, its magnetic Hamiltonian may be better ap- 
proximated by an alternating chain than a ladder, with 
the alternating chain directed along the rung of the lad- 
der (c/. Johnston et al. |Q). Our results on VODPO 
confirm the importance of PO4 pathways, and our recent 
study of magnetic excitations in (VO)2P207 single crys- 
tals [0 appears to support an alternating-chain model. 

In summary, we have used inelastic neutron scatter- 
ing on a powder sample to study magnetic excitations 
in VODPO4 • I/2D2O. Our results support a magnetic 
dimer model for this system, with an exchange constant 
of J = 7.81(4) meV, close to the value found in a sus- 
ceptibility fit. The V-V spacing within the dimer how- 
ever is found to be 4.43(7) A, which indicates that the 
dimer pair has been misidentified in previous work. The 
dimer pair we identify has a V-O-P-O-V superexchange 
pathway through a PO4 group, which evidently leads to 
a stronger magnetic interaction than the shorter V-O-V 
pathway. 

Strong superexchange through covalently bonded inor- 
ganic complexes is presumably important in many mag- 
netic insulators. For example, an exchange of k, 30 meV 
through WOg complexes has recently been identified in 
CUWO4 [|ll| and -0-Mo-O- pathways have been noted 
in FeMo04Cl jl^. The various topologies of covalent 
groups should allow the synthesis of a wide range of low- 
dimensional magnetic materials 
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Fig.l. (color) (a) The structure of VODPO at lOK 
is shown. Edge-sharing VOe octahedra (green) are con- 
nected by PO4 tetrahedra (red). These layers are sepa- 
rated by D atoms (grey). An orthorhombic unit cell is 
indicated by the dashed Hne, and the crystal axes are in- 
dicated at lower left, (b) This shows a schematic projec- 
tion of the crystal structure in the a — b plane illustrating 
possible superexchange pathways between ions. V 
ions are shown in green, P in red, and O in blue. PI is 
a V-O-V pathway (dvv=3.09 A), and P2 is V-O-P-O-V 
(dvv=4.32 A). 

Fig. 2. Scans in energy at a constant wavevector trans- 
fer of Q = 1.0 A~^ are displayed. The data are scaled to 
a fixed number of incident neutrons, as determined by a 
monitor counter. Monitor 100 corresponds to roughly 1 
minute of counting time per point. The error bars rep- 
resent one standard deviation. The profiles at different 
temperatures have been offset by a constant number of 
counts for clarity of presentation. The solid lines are fits 
to a Gaussian plus a sloping background as discussed in 
the text. 



Fig. 3. A scan in wavevector transfer Q at a constant 
energy transfer of 7.75 meV (essentially the dimer exci- 
tation energy) is shown, with some representative error 
bars. The solid line is a fit to the dimer model, as de- 
scribed in the text. 



Fig. 4. (color) This shows a two-dimensional color plot 
of the combined constant-Q scans. The contour lines are 
the results of the global fit described in the text. Counts 
per monitor 100 are given by the color bar. 
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